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Stationary and uniform entanglement distribution 
in qubit networks with quasi-local dissipation 
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We consider qubit networks where adjacent qubits besides interacting via XF-coupling, also dis- 
sipate into the same environment. The steady states are computed exactly for all network sizes and 
topologies, showing that they are always symmetric under permutation of network sites, leading to 
a uniform distribution of the stationary entanglement across the network. The maximum entangle- 
ment between two arbitrary qubits is shown to depend only on the total number of qubits in the 
network, and scales linearly with it. A possible physical realization by means of an array of doped 
cavities is discussed for the case of a linear chain. 

PACS numbers: 03.67.Bg, 03.65.Yz 



I. INTRODUCTION 

It has been understood for a long time that entangle- 
ment represents a quintessential and characteristic trait 
of quantum mechanics [l|. Quantum phenomena are by 
now very well known to be key resources for communica- 
tion and computation [2], and it has been recently ques- 
tioned whether theyplay a functional role in certain bio- 
logical processes [j, |j| • Due to its fragility under environ- 
ment induced decoherence, entanglement is commonly 
considered to be an elusive physical phenomenon that 
can observed only in the most elementary systems and 
on the shortest time scales. Nevertheless, together with 
a variety of entanglement preserving mechanisms that 
have been put forward [5l-KL3|. the idea is now spreading 
that it can persist on relatively long time scales, even in a 
noisy environment, if suitable conditions are fulfilled. To 
achieve stationary entanglement in spin systems it is suf- 
ficient to have quasi-local (two-body interaction) Hamil- 
tonian and local dissipation [7[ , or local Hamiltonian and 
quasi- local (two-body) dissipation |8l-[l0|. Till now, these 
two possibilities have been studied separately or for sys- 
tems composed of a small number of qubits. The main 
aim of this paper is to study the effects of both quasi- 
local interaction and dissipation in a system composed 
of an arbitrary number of qubits. Our goal is to deter- 
mine general conditions for stationary entanglement and 
characterize its distribution among qubits. 

We consider a family of models of quantum net- 
works consisting of n qubits with onsite energy and 
XY interaction between adjacent qubits. Moreover, a 
non-Hamiltonian dynamical term is added within the 
quantum master equations formalism |13| . The lat- 
ter describes quasi-local dissipation coupling of adjacent 
qubits, which can be understood as arising from the co- 



herent damping to the same, zero-temperature, bosonic 
environment. We compute the steady states for any size 
and network topology. This allows us to characterize 
the given model of quasi-local dissipation as a means for 
distributing stationary entanglement over a generic net- 
work. We found that the steady states are largely inde- 
pendent on the dynamical features of the model, like the 
strength of the onsite energy or of the XY interaction, 
and on the network topology. The steady states are al- 
ways symmetric under permutation of the network sites, 
yielding a uniform distribution of entanglement across 
the network. In particular the maximum attainable en- 
tanglement between any pair of qubits, measured by con- 
currence |14{ . equals 2/n ebits and is independent on the 
relative position of the two qubits. Furthermore, we in- 
vestigate the steady-state entanglement as a function of 
the initial state. Finally, for the special case in which the 
network reduces to a chain we discuss a possible physical 
realization by means of an array of doped optical cavities. 
The paper proceeds as follows. In Sec.|ll]we introduce 
the model of the qubit network; in Sec. IIIII we compute 
the steady states of the qubit network; the distribution of 
stationary entanglement across the network is discussed 
in Sec. IIVI in Sec. |V] a possible physical realization is 
introduced; Sec. lVIl is devoted to conclusions. 



II. THE NETWORK MODEL 

We consider a network of n qubits defined by a con- 
nected graph G, with vertices V{G) and edges E{G), 
where a qubit system is sitting at each vertex of the graph 
and the edges identify two body interactions between the 
qubits. Let 
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be the (symmetric) adjacency matrix of such a graph. 
Then, by considering on site energy and XY- 



interaction, the network Hamiltonian is defined as 

n -. n 



k=l 



k^l=l 



ffe 



(2) 

where w^ and a]., a^ denote respectively the energy and 
the raising, lowering operators of the fcth qubit. Further- 
more Jk^i is the coupling strength between qubits k and 
I. 

We assume that the dynamics of the network is de- 
scribed by the master equation {h = 1) 



where 



P = 'C(p) , 



C{p) = -i[H,p\+V{p)^ 



(3) 



(4) 



Here, adjacent qubits interact both directly through H 
and indirectly through a non-Hamiltonian term given in 
the Gorini-Kossakowski-Sudarshan-Lindblad form [13| . 



V{p)= X [A{G)]k,iVk,{p) : 

k^l=l 



(5) 



where 



T^k.iip) = ^ {^Lk.ipLli - LliLk,ip - pLliLkj ] , 



with 7fe.; > and 



Lk.i = o-fc -I- fj/ . 



(6) 



(7) 



The non-Hamiltonian term 'D{p) describes a Markovian 
damping process in which pairs of adjacent qubits coher- 
ently decay into the same zero-temperature bosonic bath, 
with decay rates ^k,i- 



III. THE NETWORK STEADY STATES 



A pure state \'ili) G C ®" is a steady state of the net- 
work if it satisfies £(|V')('0l) = 0. A characterization of 
the pure steady states of an open quantum system un- 
dergoing a Markovian dynamics jl3[ has been provided 
in [il| . Following [11[ , a pure steady state of our network 
model is characterized by the conditions: 

1. [A{G)]k,iLkM) = Afe.dV), for all k, I, with A^,, € 
C; 



2. 



iH+\Y.k^i[A{G)]k,ak,iLl,Lk,]m = A|^.), 
with A e C: 



3. Rc(A) = \Y.k^i[A{G)]k,nkA\kj?, where Re(A) 
denotes the real part of A. 



To compute the pure steady states of the network we 
first notice that the operators Lk^i = Uk + <^i are nilpo- 
tcnt, hence admitting only vanishing eigenvalues. Thus, 
condition [T] reads 

[A(G)]fcjK+aOIV^)-0. (8) 

A pure state can be expanded in the standard basis, 

IV') = X! V'ai,...,a„|ai,...,a„), (9) 

ai ,...,a„— 0,1 

where (7^.crfc|ai, . . . , a„) = flfejai, . . . , a„). For any pair of 
adjacent sites, k, I, Eq. ([5|) implies 

0= 2J '0ai,...U,...,a,,...,a„|ai,. . . ,Ofc, . . . ,a/,. . . ,a„) 
ak.ai=Q,l 

+ V'ai....,afc,....li,....a„|ai, ■ • ■ , Ofc , • ■ • , 0; , . . . , a„) , 

(10) 

where the notations 0^, Ik, 0/, 1/ are used to indicate 
that afc = 0, Ofc = 1, a/ = 0, a/ = 1, respectively. This in 
turn yields 



V'ai,...,lfc,---,0,,...,a„ + ■0ai,...,Ofc 
V'ai,...lfc, ...,!,, ...,«,. =0- 



0, (11) 

(12) 



If there arc no isolated points in the network, these con- 
ditions imply that the pure steady states can contain 
at most one excitation. They can be written as super- 
position of the network vacuum, |0) = |0i, . . . , 0„), and 
the single excitation states, |fc) = |0i, . . . , 1^, . . . , 0„) for 
fc = 1, . . . , n. To simplify the notation we can expand the 
pure steady states as 



= a|0)+/3X^fe|fc), 



fc=i 



where the condition (|lll) implies 

[A{G)]k,i i^k + i'l) = ■ 



(13) 



(14) 



We can now distinguish two situations according to the 
network topology: 

(i) If the network does not contain cycles or it contains 
only cycles with a even number of edges, the solu- 
tion is given by ipk — (~l)"'"0ij where Uk is the 
number of edges connecting the fcth site with the 
first one; 

(ii) Otherwise, if the network contains cycles with odd 
number of edges, the only solution is obtained by 
putting ■0fc = for all k. 

In conclusion, we get that the pure steady states have 
the form 



|V)=a|0)+/3|H) 



(15) 



where 



|H) 



1 " 



(16) 



k=l 



The coefficients a, {3 £ C are arbitrary if the network 
topology fulffils (i). On the other hand, we have to put 
/3 = if (ii) holds. Being interested in the distribution 
of stationary entanglement, in the following we assume 
that (i) is verified. 

Then, the condition [2] reads 



iiI(a|0)+/3|H)) = A(a|0)+/3|H)) , 
which may have two independent solutions: 



(17) 



• The first solution is obtained for a = under the 
conditions 



iX = ujk-J2[MG)]k,iJk,i, 



for all fc = 1, . . . , n. 

• The second solution is obtained for /3 
A = 0: 



(18) 



0, with 



In the degenerate case, uJk — J2i[-^{G)]k,iJk,i = 0, 
there exists a two-dimensional steady subspace Hg = 
span{|0),|H)}. Otherwise if ujk - Ei[MG)hiJk,i = 
const. 7^ 0, the only pure steady states are |0) and |H). 
Finally, we notice that the condition [3] is satisfied in both 
cases, since Re(A) = 0. 

Furthermore, it is worth noticing that it could be pos- 
sible that other mixed steady states exist, which are not 
in the convex hull of pure steady states. 

The steady states of our models do fulfill 
[A{G)]k,iLk,i\ilj) = 0. The steady states satisfying 
such a property are called dark states. A uniqueness 
theorem for the dark states has been provided in [ll|, 
for a system admitting a subspace of dark states. We 
can hence apply this theorem in the degenerate case, in 
which the subspace Hs = span{|0), |H)} is a subspace of 
dark states. According to this result, if there exists no 
subspace S with S J- Hs such that [A{G)]k,iLkjS C S, 
then the only mixed steady-states are in the convex hull 
of Hs- It is easy to show that a subspace with such a 
property does not exist for our models. To show that it 
is sufficient to notice that since the operators L^^i are 
nilpotent, a subspace S which is invariant under the 
action of all the L^j's must necessarily include |0) or 
|H). Thus, S cannot be orthogonal to Hs- 



IV. STEADY-STATE ENTANGLEMENT 
DISTRIBUTION 

In this section we study the steady-state entanglement 
between two arbitrary network qubits. Due to the sym- 
metric form of (|16p , the extension of the analysis of mul- 
tiqubit entanglement is straightforward. 



Since we are interested in the distribution of steady- 
state entanglement across the network, in the following 
we will restrict to the case in which both condition (ii) 
and Eq. (fT5|) are fulfilled. Under this hypothesis, the 
reduced steady state of two arbitrary qubits at site k 
and j of the network is necessarily of the form 



n/^J - f 1 _ ^ ) |00)(00| + ^|<J'^^J)(*'^^J| 



(19) 



where \^'''^} = [|10) + (-1)"j-"'=|01)] /v^ is a maxi- 
mally entangled state. It is worth noticing that the re- 
duced steady-state contains only one free parameter, p. 
Such a parameter is determined by the initial state of the 
network. To fix the ideas, we consider the concurrence 
jl4| as an entanglement measure for the reduced state of 
the two qubits. The stationary concurrence of the two- 
qubit reduced state is as well a function of p and the total 
number of qubits, namely 



C =^ 



(20) 



In order to evaluate p for a given initial state of the 
network, let us notice that 

n 

mmm = 2 E [A{G)]k,ik.mL,Mt)LUm , 

(21) 
where we have used the fact that [A(G)]fejLfc_i|H) = 0, 
and that the vectors Lj, JH) are superpositions of states 
containing two or more excitations. Moreover, we re- 
mark that the total number of excitations in the net- 
work cannot increase under the evolution dictated by 
the master equation ((4]). Therefore, we conclude that 
if the initial state p{0) contains up to one excitation, 
it follows that (H|Lfe,ip(t)L^ JH) = for any t > 0, 
which in turn yields that the quantity (H|p(i)|H) is a 
constant of motion. Eventually we get p = (H|/9(0)|H), 
which allows us to compute the steady-state parameter 
p for any initial state of the network, provided it con- 
tains up to one excitation. Let us further explore this 
setting by assuming that the network is initialized in a 
state containing a single excitation over m qubits, that 
is, IV'(O)) = J2T=i '^jl^j)- The maximum stationary con- 
currence of the two-qubit reduced state is hence obtained 
by maximizing p = KHj-f/'o)!^- It follows that the optimal 
choice for the initial state is 



1^(0)) 






(22) 



j=i 



yielding p = m/n. 

We define Cs(1,to) as the maximum stationary con- 
currence that can be achieved by preparing the network 
into an initial state containing up to 1 excitation over m 
qubits. We then have obtained that 



CJl,m) 



2m 



(23) 



To go beyond the single-excitation setting, we have ana- 
lyzed numerically the achievable stationary concurrence 
for initial states containing more than one excitation. By 
defining Cg {N, m) as the maximum stationary concur- 
rence for an initial state containing up to N excitations 
over m qubits, our numerical investigations suggest to 
conjecture that 



(24) 



a(7V,m)<C,(l,m)^ 



where the optimal network initial state is as in Eq. ([22]) . 
that is, a single-excitation initial state is sufficient to 
achieve the overall maximum concurrence for a given m. 



V. PHYSICAL REALIZATION 

We sketch here a possible physical realization by an 
array of n doped cavities coupled via optical fibers (see, 
e-g-, [la llal)- Actually, we restrict our attention to the 
case in which the network is a linear chain with open 
boundary conditions. The case of periodic boundary con- 
ditions can be analyzed in a similar way. In the case of 
a linear chain, the model generalizes that introduced in 
[lO| where the Hamiltonian term is dropped. 

Each cavity is doped with a two-level atom and is cou- 
pled by optical fibers to the next-nearest cavities. We 
denote as Cfc, c^. the ladder operators of the A:th cavity, 
coupled to the levels \g)k, \e)k of the fcth atom. Neigh- 
boring cavities are in turn coupled through a single fiber 
mode, having ladder operators a^, a^. Furthermore, 
we assume that the kth fiber mode interacts with its 
bosonic environment, described by a collection of oper- 
ators {bk,j,bl, }. The Hamiltonian of the system in the 
rotating wave approximation is given by 



H = Hf'""' + W 



(25) 



where 



Hf^^^ ^Y.^lc,ck 



J , 



^ujl\e)k{e\ + ^a;{4afc 



k=l 



k=l 

n-1 

fc=i j 



fc=i 



(26) 



and 



7L 

n-l 

y^ Jk Wk (cfc + Cfe+i) + H.c. 



(27) 



fc=i 

n-l 

+ ^^?7fc,j [albk,j+H.c. 
k=i J 

The first and second term in if /'''^^ are the free Hamil- 
tonian of the cavity field and the two level atom inside 



each cavity, the third and forth term describe the free 
Hamiltonian of the fibers modes and of the environment 



of each fibers with mode frequencies uj\ 



k\'^k^ 



J 



and 



respectively. Also, the first term in the i/™* describes the 
interaction between the cavity mode and the atom inside 
the cavity with the coupling strength Jk , the second term 
is the interaction between the cavity and the fibers modes 
with the coupling strength Jk and the last term is the in- 
teraction between the fibers and their bosonic baths with 
the coupling strength rjk^j. To write the above Hamilto- 
nian we assumed the cavities are in the strong coupling 
regime, i.e., / 3> k", k'^, where k° and k'^ are the atomic 
and cavity decay rates respectively. So, we assume that 
both the decay rates are negligible compared with the 
coupling between the fibers and their environments. 

The first two terms of H^'^^'^ and the first term of 
iJ*"* can be jointly diagonalized in the basis of polari- 
tons [la]. On the resonance between atom and cavity, 
i.e., uj'f. = ojI = ujk, the polaritonic states \n±)k = 
(|g,n)fc±|e,n-l)fe)/V2, with energies -B^ = nujk±fkV^, 
are "created" by the operators P^" ' ~ |7i±)a;(,9, 0|. 
Due to photon blockade and in the Mott phase, dou- 
ble or higher occupancy of the polaritonic states is 
prohibited, hence the only states to be considered are 
|l,±)fc, with energies Wfc ± fk [13, Um- Moreover, 
in rotating-wave approximation and interaction picture 
the inter-converting terms between different polaritons 



P, 



(l-)tp(H-) 



5(l+)tp(l-)^ 



f.,^ and Pj: -Pfc+i ) in the interaction Hamil- 



tonian are fast rotating and they average to zero. So, if 

initially polaritons are created solely by Pp. , which is 
possible by applying a global external laser to the atom- 
cavity system [16], the polaritonic state ]1, -t-) will never 
appear. Then Hamiltonian ()25p . taking into account that 
each polariton can be treated as a two level system with 



ladder operator a. 



11, 



)fc(g, 0], can be rewritten as 



ii = ^(wfc - /fc)4^fe + E '^fe^Ofc + E E '^lAj^kj 

k=l fc=l fe=l j 

fe=l 

fe=i j 

By adiabatic elimination of the fibers mode operators 
we obtain the effective Hamiltonian 



h'"=e-^4-^+j:j:<a^ 



'k,j 



fe=l 
n-1 



fc=l j 



^ J'k ((^Uk+i + <yi+i<yk 



(29) 



fe=i 

n-l 



E E ^^J [^iM + ^'^-+1) + ^■' 
fc=l J 



with 



2^1 



J2 J2 

^'k—^k-fk J- H jSk.l H j5k,n , 






72 



2J.^ 



%,j 



Jk1]k,j 



Wi 



(30) 
(31) 
(32) 
(33) 



This Hamiltonian describes a qubit chain with XY inter- 
action where nearest-neighbor qubits dissipate into the 
same bosonic bath. By tracing out the bosonic baths, 
which arc assumed at zcro-tcmperaturc, and in the Born 
and Markov approximations, one finahy obtains the mas- 
ter equation dU describing the polariton system. 



decay rates). The maximal amount of steady-state en- 
tanglement that can be achieved between two arbitrary 
qubits only depends on the size of the network and de- 
creases linearly with it. The steady-state entanglement is 
also a function of the network initial state. Furthermore, 
our analytical results, supported by numerical evidences 
leads us to conjecture that the optimal initial state of 
the network is a symmetric, Dickc-likc, state containing 
a single excitation. 

An array of doped optical cavities coupled by optical 
fibers is also discussed as a physical implementation, at 
least for the case of a network reducing to a linear chain. 
Another system could be that of planar arrays of trapped 
electrons used for quantum information processing 19j . 

Finally, the performed study lends itself to consider 
extension from 2-body dissipation to n^-body dissipation 
in n qubits network and to analyze the scaling properties 
of entanglement vs rid/n. This is left for future explo- 
rations. 



VI. CONCLUSION 

We have presented a characterization of the steady- 
states of qubit networks where adjacent qubits are cou- 
pled both directly via a XY interaction, and indirectly 
via the coherent dissipation into the same bosonic bath at 
zero temperature. We have determined conditions allow- 
ing the distribution of steady-state entanglement. Rather 
interesting, the features of the steady-state entanglement 
are largely independent on the network topology and 
on the dynamical details (e.g., coupling constants and 
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